La NMR investigation in underdoped Lai .93S10.07CUO4 
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We report 139 La and 63 Cu nuclear magnetic and quadrupole resonance (NMR/NQR) studies in 
an underdoped Lai.93Sro.o7Cu04 single crystal, focusing on the 139 La NMR in the normal state. 
We demonstrate that the local structural distortions in the low-temperature orthorhombic structure 
cause the tilting of the direction of the electric field gradient (EFG) at the nuclei from the c axis, re- 
sulting in two NMR central transition spectra at both the 139 La and 63 Cu nuclei in an external field. 
Taking into account the tilt angle of the EFG, the temperature dependence of the 139 La spectra 
allowed us to determine the 139 La Knight shift and the structural order parameter. The angle and 
temperature dependence of the 139 La spectrum is in perfect agreement with the macroscopic average 
structure and proves a displacive transition. The 139 La nuclear spin-lattice relaxation rates, T^ 1 , 
suggest that Lai.gaSro.oyCuCU undergoes a gradual change to a temperature-independent paramag- 
netic regime in the high temperature region. Both the spectra and T7" 1 of the 139 La as a function 
of temperature reveal a sharp anomaly around Ts = 387(1) K, implying a first-order-like structural 
transition, and a dramatic change below ~ 70 K arising from collective glassy spin freezing. 



I. INTRODUCTION 

La2- 2: Sr 2: Cu04 (LSCO) is a single-layered cuprate 
where the CuC>2 planes are separated by a La2_ 2; Sr 2 ;02 
block layer. A mismatch between the Cu02 planes and 
the intcrlaycr causes a structural transition from a high 
temperature tetragonal (HTT) to a low-temperature or- 
thorhombic (LTO) phased where the CuC>6 octahcdra 
are tilted from the c axis resulting in a buckling pat- 
tern of the CuC>2 plane. Further partial substitution of a 
rare earth element (Eu or Nd) at the La site, or replace- 
ment of Sr by Ba induces another structural transition 
in La2- a: Sr a ;Cu04 from the LTO to a low-temperature 
tetragonal (LTT) phase. In the LTT structure, the local 
tilt of the CuOe octahedra persists but its direction is ro- 
tated by 45° and alternates by 90° in neighboring planes 
yielding the macroscopic tetragonal symmetry. Such a 
local structural distortion as a function of doping has 
a large influence on the magnetic and superconducting 
phases^ - — In particular, the LTT phase seems to be a 
prerequisite for static stripe order which may compete 
with superconductivity in cuprates.— ~— 

In view of the intimate coupling of the respective struc- 
ture to superconductivity, a detailed understanding of 
the local structure and its doping- and temperature- 
dependence is crucial to understand the high-T c mech- 
anism. Yet, there have been unsettled issues regarding 
the local structure such as whether the macroscopic aver- 
age structure fully corresponds to the local one involving 



a displacive structural transition 
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or to a superposi- 



tion of coherent local LTT variants involving an order- 
disorder transition *£r— 

While nuclear magnetic resonance (NMR) as a local 
probe is suitable to investigate the local structure, it is 
also extremely sensitive to inhomogeneities in the sample 
which often complicate the NMR spectra, leading to an 
erroneous analysis of the results. For example, the mag- 



netic phase separation scenario^ 7 - in LSCO with x = 0.06 
based on the splitting of the 139 La NMR spectrum turned 
out to be a misinterpretation,— although its true origin 
has not been thoroughly understood. 

In this paper, the local structure in a high qual- 
ity Lai.g3Sro.o7Cu04 single crystal is investigated via 
139 La NMR. The La site is most strongly affected by 
the local distortions associated with the structural phase 
transition , 14 ! 19 making the 139 La a suitable local probe 
for the local structure. Indeed, our 139 La NMR results 
show that the local tilting pattern of the CuOe octahc- 
dra in the LTO structure is accompanied by a tilting of 
the EFG at the nuclei, resulting in the splitting of the 
otherwise single resonance line. We confirmed that the 
63 Cu spectrum splits into two in the LTO structure as 
well, but with a much smaller tilt angle than the value 
for the 139 La. A strong temperature dependence was 
found in the splitting of the 139 La spectrum, which mea- 
sures the tilt angle of the CuOg octahedra which disap- 
pears abruptly above Tg. Our NMR data thus provide 
an evidence that the macroscopic average structure pa- 
rameters are sufficient to account for the local lattice dis- 
tortions and their temperature dependence. In addition, 
the 1 La spin-lattice relaxation rates measured in a wide 
temperature range are discussed. 



II. 



SAMPLE PREPARATION AND 
EXPERIMENTAL DETAILS 



14 K) has been 



SrC0 3 , and CuO 



The crystal of Lai.93Sro.o7Cu04 (T< 
grown from 4N materials of La 2 3 , 
using the TSFZ (Traveling Solvent Floating Zone) tech 
nique. Using this technique, which is described in detail 
elsewhere,— large crystals of several centimeters length 
can be grown under accurately controllable stable con- 
ditions (flux composition, temperature and oxygen par- 
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tial pressure), which in terms are the prerequisite for 
homogeneous crystals. The sample for the NMR mea- 
surements was cut out of an single crystalline rod after 
the growth procedure along the crystallographic a, 6, and 
c-axis using the Laue x-ray method. The outer dimen- 
sions measured 1.7 mm, 2.5 mm, and 0.8 mm in a, &, and 
c-direction. 

139 La (I = 7/2) and 63 Cu (I = 3/2) NMR experiments 
were performed on the single crystal of Lai.93Sro.o7Cu04 
in the range of temperature 4.2-420 K. The sample was 
mounted on a goniometer which allows an accurate align- 
ment of the sample along the external field. In order to 
check the quadrupolc frequency vq, the 139 La NQR spec- 
trum at the Zv Q (±| «-> ±|) transition (~ 18.9 MHz) 
and the 63 Cu NQR spectrum at the vq (±| ±|) tran- 
sition (~ 34 MHz) were measured. The NMR/NQR spec- 
tra were acquired by spin-echo signals as the frequency 
was swept through the resonance line using Hahn echo 
sequence (ir/2 — r — ir) with it/2 ~ 2 [is. 

The nuclear spin- lattice relaxation rates T-J - were 
measured at the central transition of the 139 La by the 
saturation recovery method and the relaxation of the nu- 
clear magnetization after a saturating pulse was fitted to 
the following equation: 

1 _m L _j^ /T ? L/Ti 

M{oo) _ \M + 44 

(!) 

, 75 ,-i5tm , 1225 28t/Ti \ 
+ 364 + 1716 J ' 

where M is the nuclear magnetization and a a fitting 
parameter that is ideally one. 
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FIG. 1: Angle dependence of 139 La NMR spectra measured at 
room temperature and 10.7 T. Solid (dotted) lines are theo- 
retical calculations assuming 8° (—8°) tilting of the principal 
axis of the EFG from the c axis, using K, = —0.09 % and 
vq — 6.3 MHz which were determined experimentally from 
this study. The anisotropy parameter rj — 0.19(2) was ob- 
tained from the calculations. 



III. RESULTS AND DISCUSSION 
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A. 139 La NMR spectra 

Figure[T]shows the central transition NMR spectrum of 
the 139 La at 300 K and 10.7 T as a function of the angle 9 
between the crystallographic c axis and the external field 
H . The 139 La spectrum consists of two resonance lines 
which both exhibit the same angle dependence with re- 
spect to H: the peaks are centered at 9^ = ±8°, respec- 
tively, where the maximum resonance frequencies are ob- 
served. The two angle-dependent sets of the 139 La spec- 
trum which are symmetric around the c axis imply the 
existence of two different sites with the same occupation 
probability. 

Such a strong angle dependence of the central transi- 
tion is expected from second order quadrupole interaction 
of the nuclear quadrupole moment Q that is non-zero for 
the nuclear spin / > i, with the EFG. For the central 
transition (— -| O ^), while the first order quadrupole 
shift is zero, the second order quadrupole shift is given 



A^ 2nd = 



32 ln H 



(1-cos 2 9) 



21(1+1)-- 



27 



- 118/(1+1)--) cos 



(2) 



where 7„ is the nuclear gyromagnetic ratio and vq is the 
quadrupole frequency given by 2>e 2 qQ /2I(2I — \)h where 
eq is the EFG. Here, without the loss of generality, we 
ignored an anisotropy parameter rj which is non-zero in 
sites of lower than axial symmetry. From Eq. @, one 
can see that when # = 0,A^ = 0so that the reso- 
nance frequency becomes one of the maxima and directly 
yields the Knight shift since a quadrupole correction is 
not needed. 

It may be possible that the 139 La is differentiated due 
to the presence of two EFG strengths which would cause 
the different quadrupole shift at the nuclei. To check 
the possibility, we measured the 139 La NQR spectrum 
in zero field at the 3vq transition. Since a well-defined 
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FIG. 2: 139 La NQR spectrum measured at the 3uq transition 
at 40 K giving rise to vq = 6.3 MHz. The NQR spectrum from 
Julien et al.— measured in the sample of x = 0.1 is compared. 
The narrower linewidth of our NQR spectrum indicates less 
chemical inhomogeneity in the single crystal. 

single NQR line as shown in Fig. [5] indicates the unique 
value of the EFG in the system, the two 139 La peaks 
are attributed to the local tilting of the direction of the 
EFG at the 139 La with respect to the c axis^ where 9 m 
corresponds to the angle at which the direction of the 
EFG is parallel to H. 

In order to confirm this quantitatively, we performed 
exact diagonalization calculations of the nuclear Hamil- 
tonian which can be written as, 



H=- 7nS[l 



37? 



£(0)]I 

-I 2 + 



H 



'/ 



(3) 



where K is the Knight shift which is usually angle depen- 
dent and I is the nuclear spin operator. 

In the calculations, we used vq = 6.3 MHz from the 
iv Q transition (Fig.© and K = -0.09 % (Fig.©. Here 
we neglect a possible angle dependence of /C. Thus, r\ 
is the only unknown parameter in the Hamiltonian. By 
adjusting 77, we accurately reproduced the resonance fre- 
quencies of the spectra as a function of 9, drawn as solid 
(dotted) lines in Fig. [TJ From these calculations, the 
anisotropy parameter r\ = 0.19(2) was deduced. Since rj 
can absorb the effect of the angle variation of K,, the ac- 
tual value may slightly differ from the one deduced here. 
The excellent agreement of the theory with the experi- 
mental data indicates a staggered pattern of the tilting 
of the EFG at the 139 La around the c axis, due to the 
corresponding local structural distortions. 

As a matter of fact, this is well supported by the LTO 
structure in Fig. [3J which was drawn using structure pa- 
rameters for x = 0.075 at 10 K given in Ref. [l^. Clearly, 
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FIG. 3: Structure of La2-2,Sr a; Cu04 in the LTO phase at 10 
K for x — 0.075, using the structure parameters reported in 
Ref. 01, viewed along [001] (a) and along [100] (b). The unit 
cell is depicted as the dotted line. The Cu06 octahedra is 
tilted along the [010]lto (or [110]htt) direction by 4°. In 
addition, La(Sr)-La(Sr) along the c axis is tilted by 1.5° in 
opposite direction of the distortion of the CuO g octahedra, 
being consistent with a much larger tilted angle of the EFG 
at the 139 La. Note that the tilt direction of both the CuO e oc- 
tahedra and the La-La is maintained along [010]lto, whereas, 
along the diagonal direction i.e., [110]lto, it is rotated alter- 
nately by 180°, causing the so called buckling of the CUO2 
plane. 



the main axis of the CuOe octahedra is tilted by 4° along 
the [010] direction (or [110] in the HTT setting). In ad- 
dition, the position of La(Sr) with respect to the c axis 
is slightly shifted in opposite direction of the distortion 
of the CuOg octahedra, which may cause a larger tilting 
of the EFG at the 139 La than at the 63 Cu. Furthermore, 
the fact that the tilt direction of the octahedra is reversed 
alternately along the [110]lto, which induces the buck- 
ling of the Cu02 planes^ accounts for the two 139 La 
peaks whose intensities and positions as a function of 9 
are symmetric around the c axis. 

If the two 139 La NMR lines arc the result of the dis- 
torted CuOg octahedra in the LTO phase, one should 
expect that there are two Cu sites which are distinguish- 
able in field as well. However, the presence of the two 
central transitions of the 63 Cu has not been reported so 
far. It may be possible that the tilt angle of the EFG 
at the 63 Cu is too small to be detected, especially with 
the relatively larger linewidth of the 63 Cu NMR spec- 
trum than that of the 139 La. In order to check whether 
there exist two 63 Cu NMR lines in the LTO phase, the 
63 Cu spectrum has been carefully examined by varying 
9 at 300 K. Indeed, as shown in Fig. the angle depen- 
dence of the 63 Cu spectrum unravels two central transi- 
tions which are centered at 9^ = ±2° respectively, being 
symmetric around the c axis. The solid (dotted) lines 
are from theoretical calculations with Eq. (© using the 
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FIG. 4: Angle dependence of the 63 Cu NMR spectra measured 
at 300 K and 8.2 T. The solid (dotted) lines are theoretical 
calculations assuming 2° (—2°) tilting of the direction of the 
EFG at the 63 Cu from the c axis, with parameters vq = 34 
MHz, K = 1.18 %, and r\ = 0. 



values of parameters, vq = 34 MHz and /C c = 1.18 %. 
Unlike the case of the 139 La, a non-zero 77 is not needed 
to fit the data. The negligible 77 and the much smaller 
6 m than that of the 139 La are also consistent with the 
LTO structure depicted in Fig. [31 where the La itself is 
displaced. Therefore, we conclude that the two peaks at 
both the 139 La and the 63 Cu are a natural consequence 
of the HTT-LTO structural transition in La2_ a: Sr a ;Cu04. 

We are aware that the qualitatively similar argument 
was given by Julicn et alJ£ where they observed four 
139 La peaks which were attributed to the tilting of the 
direction of the EFG about the c axis. However, the 
four peaks with various signal intensities and linewidths 
are difficult to understand in detail and are not in ac- 
cordance with only two different La sites as expected for 
the LTO structure. We conjecture that the sample used 
by Julien et al. consisted of two or more domains which 
have slightly different orientations. 
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FIG. 5: Temperature-dependence of the 139 La NMR spectra 
in an external held tilted by 9 ~ 8° from the crystallographic c 
axis. In this case, pi corresponds to the line in which the EFG 
at the 139 La is aligned along H, giving rise to the Knight shift 
for H || c since there is no quadrupole shift. At 387 K, the 
two peaks collapse to a single peak, indicating the LTO-HTT 
structural transition. At low temperatures, 139 La spectrum 
is significantly broadened, implying the slowing down of the 
spin fluctuations due to the spin-density wave (SDW) order. 
Each spectrum was scaled with temperature to remove the 
Boltzmann contribution to the intensity. 



B. La Knight shift and structural order 
parameter 

Next we discuss the temperature dependence of the 
139 La spectra. In order to extract the Knight shift 
for / > 1/2 in non-axial symmetry, the second order 
quadrupole effect which shifts the central transition, as 
given by Eq. (J2J, should be taken into careful consider- 
ation, particularly for a very small Knight shift which is 
the case for the 139 La NMR in La 2 - a; Sr K Cu04. Taking 
advantage of the fact that the second order quadrupole 
shift vanishes when H is applied along the direction of 
the EFG, we deliberately rotated the sample with respect 
to the c axis by 8°. Then, the second order quadrupole 
shift for the peak pi is zero and the intrinsic Knight shift 
function of temperature could be extracted. 

With increasing temperature from 300 K, the sepa- 
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ration between the two peaks becomes smaller and the 
spectra collapse abruptly into single spectrum at Tg = 
387(1) K, the structural transition temperature Tg from 
the HTT to the LTO phased This is another direct 
evidence that the two 139 La central peaks arise from the 
tilted CuC>6 octahedra associated with the LTO struc- 
ture. The T-dcpendcncc of the resonance frequencies of 
the two lines is shown in Fig. [BJa) . The single line above 
Tg is split into two at Tg, in which one line (pi) shifts 
up and the other (p2) shifts down in frequency from the 
single line above Tg. The line pi gives rise to the Knight 
shift for H | c, JC C , which is drawn in Fig.^b). Note that 
the Knight shift data above Tg is obtained after rotating 
back the sample along the c axis, since the EFG should 
be directed along the crystallographic c axis in the HTT 
structure due to the axial symmetry with respect to the 
c axis. 

Additional valuable information was obtained from 
Fig-EIa). Just below Tg, the separation between pi and 
p2 rapidly increases but approaches a constant value at 
low temperatures, that is, the separation between the 
two lines Av as a function of temperature behaves as the 
structural order parameter as shown in Fig.[BJb). At Tg, 
the transition is clearly discontinuous, indicating a first 
order nature of the transition, in contrast to second-order 
character reported in the diffraction measurement^ On 
the other hand, except for the discontinuous change at 
the transition, the temperature dependence of Av shows 
the behavior expected in a second order transition. 

The narrow single resonance line in the HTT phase in- 
dicates that within the resolution of the NMR experiment 
the local tilting of the CuOg octahedra vanishes immedi- 
ately above Tg rendering the local tetragonal symmetry, 
in contrast to the experimental reports that the LTO- 
type local tilt remains even in the HTT phase d 4 ' 24 Fur- 
thermore, the abrupt splitting of the spectrum at Tg con- 
firms a displacive transition, not an order-disorder one. 

Since Av is a measure of the tilt angle of the EFG at 
the nuclei, its rapid suppression with increasing temper- 
ature from 300 K toward Tg may cause a finite second 
order quadrupolc shift affecting the 139 La Knight shift 
which was obtained with the tilt angle 8° determined at 
300 K. Then, the actual Knight shift data just below Tg 
may be slightly enhanced, in such a way that data are 
smoothly connected to those in the HTT phase. Like- 
wise, Av itself should be affected with a small error. 



C. La spin-lattice relaxation rate 

Figure [7] shows the spin lattice relaxation rate T± 1 of 
the 139 La as a function of temperature in the range 4.2- 
420 K. With increasing temperature from room temper- 
ature, T^ 1 is rapidly enhanced. We found a discontinu- 
ous jump at Tg = 387 K which is followed immediately 
by a steep drop above Tg. The sharp anomaly of T^ 1 
around Tg is consistent with the drastic change of the 
139 La spectrum at Tg (Fig. [5]) and supports a first order 
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FIG. 6: (a) Temperature-dependence of resonance frequencies 
of two NMR peaks in Fig. [5] (b) Knight shift of the 139 La is 
extracted from the pi line. The separation of the two peaks 
Av(T) behaves as a structural order parameter. Av reveals 
first order-like discontinuous transition at Ts = 387 K, while 
it exhibits the second order-like temperature dependence ex- 
cept the region of the transition. Solid line is a guide to eyes. 
Below 70 K, data deviate due to the onset of the magnetic 
broadening associated with the incommensurate SDW order. 



structural transition. Upon further increasing the tem- 
perature, T-j - continues to increase up to 420 K which is 
a limitation of our equipment. 

Interestingly, the temperature dependence of ^J 1 " 1 
resembles that observed in Lais-zEuo^Sr^CuO^i^ 
where 139 Tf 1 rapidly drops below the LTT-LTO struc- 
tural transition. The HTT-LTO transition, however, is 
clearly manifested through an anomaly in the narrow 
range of temperature which is simply added on top of 
the rapid drop of Tj~ . Note that T-f already started to 
decrease at a much higher temperature than Tg. 

Rather, we conjecture that the rapid upturn of T-j - 
is associated with a gradual change from a system with 
a local moment character to a more itinerant system, 

Tr 1 of 



as argued by Imai et al£± They report that ± 1 
the 63 Cu forms a broad maximum near 200 K in sim- 
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inated by the doped hole carriers in this temperature 
range. Therefore, we ascribe the rapid drop of Tr 1 at 
high temperatures to an evolution from a paramagnetic 
regime at T > J cx /2 to a regime below 300 K where the 
scattering with doped holes dominates the spin fluctua- 
tion properties, and thus the 139 La Tr . 

Below ~ 70 K, T^ 1 is enhanced in an extremely fast 
rate, which is ascribed to the critical slowing down of spin 
fluctuations associated with glassy spin freezing. Such a 
rapid upturn of T-T 1 is compatible with the significant 
broadening of the 139 La NMR spectrum in the same tem- 
perature range as shown in Fig. [5] Down to 4.2 K, a 
maximum of Tf 1 , which gives a characteristic spin freez- 
ing temperature T g , was not identified. This is somewhat 
inconsistent with previous 139 La NMR results^^ and 
suggests a reduced T g which is possibly related to less 
disorder in our sample. 



FIG. 7: 139 La T^ 1 as a function of temperature. A very 
sharp anomaly at Ts = 387 K was detected with the transi- 
tion width less than 10 K. Below 70 K, the critical slowing 
down due to the SDW order causes extremely steep upturn of 
T[~ . Inset shows (TiT) -1 of the 139 La versus temperature. 
Note that (TiT) -1 increases more than four decades as the 
temperature increases from ~ 100 K down to 4.2 K. 

ilarly Sr-doped samples, approaching a temperature in- 
dependent value 63 ~ 2500 s" 1 at sufficiently high 
temperatures! 26 ' 27 If the temperature and doping inde- 
pendent 63 at T ~ Jcx/2, where J ex is the exchange 
coupling constant, imply the exchange narrowed para- 
magnetic limit, one may expect that 139 should be 
dominated by J cx and the hyperfine coupling Am accord- 
ing to the relation ; 28 ' 29 

— oc A "^TT), (4) 

T\ oo Jcx 

where we neglect an anisotropy of Am for simplicity. A 
very rough estimation from the Knight shift values, which 
are proportional to Am, obtained at room temperature 
gives 139 T{^ = 63 T^( 139 /C/ 63 /C) 2 ~ 10 s^ 1 which is 
surprisingly close to the experimental values near 400 K. 
Then, one may anticipate that with increasing tempera- 
ture beyond 420 K i 39 ^ -1 increases further and decreases 
ultimately, or gradually saturates, to a constant value 
that is close to ~ 10 s _1 , being temperature-independent. 
This indicates that the 139 La relaxation rates at high 
temperatures near J ex /2 are governed by the same local 
spin dynamics probed in the 63 Cu results? 2 ^ Therefore, 
the S — | antifcrromagnctic 2D Hciscnberg mode l 29 i 30 
may be indeed appropriate to describe the high temper- 
ature spin dynamics in LSCO even at T r * j Jcx/^ • 

Furthermore, we find that T-T 1 decreases linearly with 
temperature below 300 K down to ~ 150 K. This is 
clearly shown by T-T 1 divided by temperature, (TiT) -1 
(see the inset of Fig. [7]) . The Korringa-likc behavior sug- 
gests that the relaxation process of the 139 La is dom- 



IV. CONCLUSION 

Through 139 La and 63 Cu NMR, we proved that the lo- 
cal tilting of the CuC>6 octahedra along [110]htt occurs 
in the LTO phase, whereas it is absent in the HTT phase 
in Lai.g3Sro.o7Cu04. The tilt direction of the CuOg oc- 
tahedra alternates along the same direction resulting in 
the buckling of the C11O2 plane. The structural order 
parameter representing the tilt angle of the CuC>6 oc- 
tahedra was successfully obtained from the splitting of 
the 139 La resonance line as a function of temperature. 
The structural transition at Ts = 387 K is weakly first 
order and displacive i.e., the order parameter is reduced 
rapidly with increasing temperature approaching Ts, and 
abruptly vanishes in the HTT phase. In other words, the 
average HTT as well as the LTO phases are fully com- 
patible with their local structures. Our data also rule out 
the possibility of an inhomogencous mixture of the two 
phases in the whole temperature range measured in this 
study. 

The 139 La spin-lattice relaxation rates, Tj - , exhibits a 
sharp anomaly at Ts, which is consistent with the abrupt 
change of the spectrum. The temperature dependence 
of T^ 1 in the high temperature region suggests that a 
temperature-independent paramagnetic regime at high 
temperatures changes gradually to a regime where the 
scattering with doped holes dominates the 139 La Tj~ . 
The strong enhancement of Tj -1 below ~ 70 K is as- 
cribed to the critical slowing down of spin fluctuations 
associated with glassy spin freezing. 

Acknowledgement 

The authors wish to thank Markus Hiicker for his help- 
ful discussions and communications. This work has been 
supported by the DFG through FOR 538 (Grant No. 
BU887/4 and No. ER342/1-3) and SPP1458 (Grant No. 
GR3330/2). 



7 



* |sbaek.fuQgmail.c om| 

1 B. Keimer, N. Belk, R. J. Birgeneau, A. Cassanho, C. Y. 
Chen, M. Greven, M. A. Kastner, A. Aharony, Y. Endoh, 
R. W. Erwin, et al., Phys. Rev. B 46, 14034 (1992). 

2 J. D. Axe and M. K. Crawford, J. Low Temp. Phys. 95, 
271 (1994). 

3 B. Biichner, M. Breuer, A. Freimuth, and A. P. Kampf, 
Phys. Rev. Lett. 73, 1841 (1994). 

4 A. Bianconi, N. L. Saini, A. Lanzara, M. Missori, T. Ros- 
setti, H. Oyanagi, H. Yamaguchi, K. Oka, and T. Ito, Phys. 
Rev. Lett. 76, 3412 (1996). 

5 B. J. Suh, P. C. Hammel, M. Hiicker, and B. Biichner, 
Phys. Rev. B 59, R3952 (1999). 

6 H.-H. Klauss, W. Wagener, M. Hillberg, W. Kopmann, 
H. Waif, F. J. Litterst, M. Hiicker, and B. Biichner, Phys. 
Rev. Lett. 85, 4590 (2000). 

7 J. M. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Nakamura, 
and S. Uchida, Nature 375, 561 (1995). 

8 H.-J. Grafe, N. J. Curro, M. Hiicker, and B. Biichner, Phys. 
Rev. Lett. 96, 017002 (2006). 

9 M. Vojta, Adv. Phys. 58, 699 (2009). 

10 C. Friedrich, B. Biichner, M. M. Abd-Elmeguid, and 
H. Micklitz, Phys. Rev. B 54, R800 (1996). 

11 M. Braden, M. Meven, W. Reichardt, L. Pintschovius, 
M. T. Fernandez-Diaz, G. Heger, F. Nakamura, and T. Fu- 
jita, Phys. Rev. B 63, 140510 (2001). 

12 B. Simovic, M. Hiicker, P. C. Hammel, B. Biichner, U. Am- 
merahl, and A. Revcolevschi, Phys. Rev. B 67, 224508 
(2003). 

13 S. J. L. Billinge, G. H. Kwei, and H. Takagi, Phys. Rev. 
Lett. 72, 2282 (1994). 

14 D. Haskel, E. A. Stern, D. G. Hinks, A. W. Mitchell, J. D. 
Jorgensen, and J. I. Budnick, Phys. Rev. Lett. 76, 439 
(1996). 

15 D. Haskel, E. A. Stern, F. Dogan, and A. R. Moodenbaugh, 
Phys. Rev. B 61, 7055 (2000). 



S.-W. Han, E. A. Stern, D. Haskel, and A. R. Mooden- 
baugh, Phys. Rev. B 66, 094101 (2002). 
M.-H. Julien, F. Borsa, P. Carretta, M. Horvatic, 
C. Berthier, and C. T. Lin, Phys. Rev. Lett. 83, 604 (1999). 
M.-H. Julien, A. Campana, A. Rigamonti, P. Carretta, 

F. Borsa, P. Kuhns, A. P. Reyes, W. G. Moulton, M. Hor- 
vatic, C. Berthier, et al., Phys. Rev. B 63, 144508 (2001). 
P. G. Radaelli, D. G. Hinks, A. W. Mitchell, B. A. Hunter, 
J. L. Wagner, B. Dabrowski, K. G. Vandervoort, H. K. 
Viswanathan, and J. D. Jorgensen, Phys. Rev. B 49, 4163 
(1994). 

M. Lambacher, T. Helm, M. Kartsovnik, and A. Erb, Eur. 
Phys. J. Special Topics 188, 61 (2010). 

G. C. Carter, L. H. Bennett, and D. J. Kahan, Metallic 
shift in NMR (Pergamon, New York, 1977). 

P. C. Hammel, A. P. Reyes, S.-W. Cheong, Z. Fisk, and 

J. E. Schirber, Phys. Rev. Lett. 71, 440 (1993). 

T. Nakano, M. Oda, C. Manabe, N. Momono, Y. Miura, 

and M. Ido, Phys. Rev. B 49, 16000 (1994). 

S. Wakimoto, H. Kimura, M. Fujita, K. Yamada, Y. Noda, 

G. Shirane, G. Gu, H. Kim, and R. J. Birgeneau, J. Phys. 

Soc. Jpn. 75, 074714 (2006). 

N. J. Curro, P. C. Hammel, B. J. Suh, M. Hiicker, 
B. Biichner, U. Ammerahl, and A. Revcolevschi, Phys. 
Rev. Lett. 85, 642 (2000). 

T. Imai, C. P. Slichter, K. Yoshimura, and K. Kosuge, 
Phys. Rev. Lett. 70, 1002 (1993). 

S. Fujiyama, Y. Itoh, H. Yasuoka, and Y. Ueda, J. Phys. 

Soc. Jpn. 66, 2864 (1997). 

T. Moriya, Prog. Theor. Phys. 16, 641 (1956). 

S. Chakravarty and R. Orbach, Phys. Rev. Lett. 64, 224 

(1990). 

M. P. Gelfand and R. R. P. Singh, Phys. Rev. B 47, 14413 
(1993). 



